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Abstract 

Using the recently developed method we calculate the crystal field parameters in yttrium and lutetium 
aluminum garnets doped with seven trivalent Kramers rare-earth ions. We then insert calculated parameters 
into the atomic-like Hamiltonian taking into account the electron-electron, spin-orbit and Zeeman interac¬ 
tions and determine the multiplet splitting by the crystal field as well as magnetic g tensors. We compare 
calculated results with available experimental data. 
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I. INTRODUCTION 


To explain certain optical and magnetic properties of rare-earth (R) materials, determination of 
crystal field parameters (CFP) is essential. When snfficient experimental data are available CFP are 
nsnally determined by the least sqnares £t. Since the nnmber of nonzero CFP depends on the site 
symmetry and may be as high as 27, such method often ends up being ambiguous. As a result, there 
has been a continuous effort to estimate CFP theoretically (for the review of various methods, see 
Ref. flfl). 


Recently, a novel theoretical approach to calculate CFP has been proposed j^. An original mo¬ 
tivation of the work was to explain the magnetic properties of rare-earth cobaltites RCoO^ where 
available experimental data do not suffice to estimate CFP. The method starts with the density 
functional theory (DFT) based band structure calculation, followed by a transformation of the Bloch 
to Wannier basis. The local Hamiltonian is then expanded in terms of the spherical tensor operators. 
Resulting CFP are inserted in an atomic-like Hamiltonian involving the crystal held, 4/—4/ correla¬ 
tion, spin-orbit coupling and Zeeman interaction. The method does not suffer from the 4/ electron 
self-interaction (the difficult problem of DFT ab initio methods). The hybridization of the 4/ states 
with other valence orbitals is taken into account via hybridization parameter A, a single parameter 
of the method. In the recent paper jj] a relatively simple way allowing to estimate this parameter 
was suggested and applied to the RiLaFs system. 

The method has been extensively tested for rare-earth doped aluminates YAlOsiR^'*' with or¬ 


thorhombic perovskite structure [31. Remarkable agreement of calculated and experimental data was 

n n 

achieved. Application to gallates RGaOs and cobaltites RC 0 O 3 [5| as well as manganites RMnOs p 
followed. Even in these cases calculations provide a fair agreement with available experimental data. 

Rare-earth doped aluminum garnets are widely used as laser materials and scintillators due to 
which extensive experimental data are available. In recent years magneto-optical properties of some 
non-Kramers ions in garnet hosts have been studied. High magneto-optical activity observed in 
some cases is of interest in microwave anmlihers and generators. Magnetooptics of Tb^'*' and Tm^"*" 


in Y3AI5O12 (YAG) is studied in Refs. ^ 
various garnets is reported in Refs. [lOj, iU 


and Ref. 
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respectively. Ma gne 


and that of Pr^+ in YAG in Ref. 


tooptics of Eu^’*' 


12 


m 


In this work we apply the new CFP method to rare-earth doped aluminum garnets. After reviewing 
the theory and computational procedure, the next section is devoted to analysis of the problem of 
determination of the parameter A in more detail compared to Ref. [J. In the following section we 
calculate crystals held parameters. We show examples of multiplet splitting by the crystal held. We 
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focus on Kramers ions and calculation of their magnetic g factors, yet more crystal-field sensitive 
quantities, and compare them with experimental data. 


II. THEORETICAL APPROACH AND COMPUTATIONAL DETAILS 


The effective Hamiltonian describing 4/ states can be written as 


H,ff=HA + Hz + H, 


CF ■ 


( 1 ) 


where Ha is the spherically symmetric, free ion atomic-like Hamiltonian (for details see Ref. Il3h . 
while Hz and Hcf are the Zeeman interaction and crystal held Hamiltonian, respectively. In the 
Wybourne notation [l^ Hcf has the form 

■ ( 2 ) 

A:=2,4,6 q=—k 

'' ih') 

where Cg is a spherical tensor operator of rank k acting on the 4/ electrons of the R ion. The 
coefficients are the crystal held parameters. Hermiticity of Hcf requires that 

Calculation of crystal held parameters consists of four steps: 

1. Standard selfconsistent band calculation with 4/ states included in the core. The results yield 
the crystal held potential, subsequently used in the next step. 

2. The 4/ as well as oxygen 2p and 2s states are treated as the valence states in a nonselfconsistent 
calculation, all other states are moved away using the orbital shift operator. Relative position 
of 4/ and oxygen states is adjusted using the hybridization parameter A (a single parameter of 
the method). 


3. The 4/ band states are transformed to Wannier basis using the wien2wannier [15[ and wan- 
nierQO Q packages. 

4. Local 4/ Hamiltonian in the Wannier basis is extracted and expanded in a series of spherical 
tensor operators. The coefficients of expansion are the crystal held parameters. 

To perform the band structure calculations in steps 1 and 2 we used the WIEN2k package 0 
with implemented augmented plane waves -|- local orbital method. For the exchange correlation 
functional we applted the genetahaed-gtadient apptoxtaation fottn Q, We used expetnnental lattice 
parameters of Y 3 AI 5 O 12 (YAG) and LU 3 AI 5 O 12 (LuAG) [19|, but the atomic positions within the 
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unit cell were optimized for each R substitution by minimizing the atomic forces. The unit cell 
in our calculations consisted of 80 atoms. The eigenvalue problem was solved in hve points of 
the irreducible Brillouin zone and the number of basis functions was ~ 7700 (corresponding to the 
parameter i?if’maa;=6.5). The calculations were non-spin polarized. The atomic radii of R (Y,Lu), 
A1 and O were 2.3, 1.7 and 1.55, respectively. 


Once the crystal field parameters were determined we used the modihed lanthanide package 


20 | 


to solve the eigenproblem for the Hamiltonian ([T]). The results provide the multiplet splitting by the 
the crystal field. From the energy dependence on external magnetic held one extracts the g tensor 
(for more details see Refs. ^,5). 


III. HYBRIDIZATION PARAMETER 


The parameter A appears due to hybridization between the rare-earth 4/ states and the valence 


Q 


In the R 


states of its ligands. Our treatment of hybridization is briehy described in Ref, 
containing orthorhombic perovskites 3,5,^ remarkable agreement between experimentally obtained 
spectroscopic data and the calculation was obtained by hxing the value of A at 0.6 Ry. With the same 
A also the magnetism was calculated and compared with experiment. However, less experimental 
data were available and the agreement, though still satisfactory, was not as good as in case of 
spectroscopy. In RiLaFs very good agreement between optical data and calculation was obtained for 
A = 0.4 Ry jj] . In the same paper A was estimated using a charge transfer energy 


A ~ Nval - 1) - N,al), 


(3) 


where N^ai) is the total energy of the ground state of the system {n^f electrons in 4/ shell of 

R ion and N^ai electrons in the valence band), while N^ai — 1) corresponds to the excited 

state in which one of the valence electrons was transferred in the 4/ electron shell. The hybridization 
parameter thus can be calculated by performing two calculations with 4/ electrons treated as the 
core states - the first one with 4/”, N^ai, the second with N^ai — 1 electron conhgurations. 

Using the above equation we calculated A for the R in question in both YAG and LuAG. The 
results, together with the data for RiLaFs and orthorhombic perovskites RMnOs, are shown in Fig. [T] 


There are several problems connected with the above method. The hrst one is connected with the 
multiplet splitting of the 4/ levels, which is not provided by the DFT calculation. In principle this 
splitting may be obtained by the atomic-like program, we are using. However, there is a problem of 
double counting of the electron-electron correlation, which would be difficult to overcome. We can 
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FIG. 1: Dependence of the hybridization parameter A on number of the 4/ electrons calculated using eq. 
[3l Dashed and dash-and-dotted lines correspond to the A values, which were adopted for calculations in 
RMnOs and i?:LaF 3 , respectively [j, l6|. 


only estimate that corresponding nncertainty of A is on the order of the crystal held splitting i.e. 
less than ~ 0.1 eV and it will change with R. Such scatter will become relatively more important in 
R:YAG and i?:LuAG compared to i^iLaFs and i^MnOa, because A calculated from (|3]) is smaller in 
garnets (cf. Fig. [T]). 

The second problem is inherent to the open core calculations. Even though the 4/ electrons are 
well localized, small part of their density leaks out of the R atomic sphere. Understandably, this 
leakage is bigger for the electron conhguration. For i?:LuAG we calculated A for three values 

of the atomic sphere radius Rmt- The result is shown in Fig. [2J 

Finally, note that eq. ([3]) is based on the hrst order perturbation theory j^, thus higher order 
correction is needed if the hybridization is important. 


IV. RESULTS 

A. CFP and multiplet splitting 

Yttrium (Lutetium) aluminum garnets have a cubic structure belonging to the laSd space group. 
The unit cell contains eight molecular units Y 3 AI 5 O 12 (YAG) or Lu 3 Al 50 i 2 (LuAG). The R impurities 
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number of 4f electrons 


FIG. 2: ii:LuAG. Hybridization parameter calculated for three values of the atomic sphere radius Rmt- 


are located on dodecahedrally coordinated Y (Lu) sites of D 2 symmetry, and nine real parameters 
are necessary to characterize the crystal held. For the dodecahedral sites, there are six possible 
coordinate system orientations, which result in six different (but spectroscopically indistinguishable) 
crystal held parametrizations. The coordinate system to which our CFP and g tensors are referred 
is shown in Fig. El 

The magnitude of parameter A entering the calculation in the second step was varied between 
2.7-10.9 eV (0.2-0.8 Ry) with the step 1.36 eV (0.1 Ry). As A decreases, the 4/ levels get closer to 


the valence band, for A = 2.7 eV the calculation becomes 


ess stable for lighter R and it crashes for 
6 | and LaFs jd] for hxed A the CFP change 


Er and Yb. Similarly as in orthorhombic perovskites fly, 3 
smoothly with the number of 4/ electrons. In Fig. Slthis is documented for A=5.4 eV and YAG. 
As a function of the hybridization parameter, CFP also change smoothly. An example for Sm;YAG 
is displayed in Fig. [5l The CFP for both YAG and LuAG and all seven R are collected in Table 
IIV A1 They were calculated taking A=5.4 eV, which is the lowest value for which the calculation 
runs smoothly for all R. 


With the knowledge of CFP, modihed program ’lanthanide’ 20|] was used to calculate the energy 


levels. Parameters of the free ion Hamiltonian were taken from Ref. 


13 


The agreement between 

theory and experiment is generally fairly good. Three examples, namely Nd^"*", Sm^"*" and Er^"*" in 
YAG are shown in Figs. El [71 [HI For lighter R the position of energy levels only slightly depends on 
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FIG. 3: YAG structure with the highlighted dodecahedrally coordinated site. Axes a and b are parallel to 
[lIO] and [110] directions, axis c is parallel to [001] direction. 



FIG. 4: YAG. CFP as a function of the number of 4/ electrons for hybridization parameter A=5.4 eV. 
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FIG. 5; Sm:YAG. CFP as a function of the hybridization parameter A. 


the hybridization parameter, for heavier R the dependence becomes stronger 


fl. 


B. g tensors 


Axes a, b, c of the system are principal axes of g. In the experiment components ga, gb, gc are usually 
determined by analysis of the dependence of EPR spectra on the direction of external magnetic held. 
There are two sites Ri, R 2 with the c axis running along the [001] direction. Axes ai,bi, are parallel 
to [ 110 ] and [ 110 ], respectively. Axes 02,62 are obtained from ai,^by a 7 r /2 rotation around c. 
This leads to an ambiguity in assessment of ga, gt to Ri, R 2 sites [ 2 ^. On the other hand, in the 
calculation the assessment is unambiguous and it refers to the R site shown in Fig. [31 

To determine g the effective Hamiltonian (1) was diagonalized with gradually increasing external 
magnetic held B. Resulting dependence of energies was then expanded up to the second power of B. 
The linear term provides g. Comparing to energies in the zero-held, g is much more susceptible to 
parameters of the calculation, in particular to A. In Figs. 9-15 the dependence ^(A) is shown for all 
R in YAG and LuAG. As mentioned above, for small A the calculations become unreliable, which 
results in huctuations in ^(A) dependence. The calculation for A = 5.4 eV and i?:YAG, R;LuAG is 
compared with the experiment in Tables jlT] and IIIIl 















TABLE I: Crystal field parameters (in meV) in R:YAG and ii:LuAG calculated for hybridization parameter 
A=5.4 eV. 


k q 

Ce:YAG 

Nd:YAG 

Sm:YAG 

Gd:YAG 

Dy:YAG 

Er:YAG 

Yb:YAG 

2 0 

-21.2 

32.8 

37.0 

52.4 

67.4 

113.1 

72.2 

2 2 

-46.6 

-28.7 

-23.8 

-14.4 

-5.6 

15.4 

2.1 

4 0 

-5.6 

-5.8 

-1.1 

-5.8 

-8.6 

-22.9 

-3.0 

4 2 

313.8 

253.9 

223.2 

206.3 

193.1 

183.8 

163.6 

4 4 

-203.9 

-161.9 

-139.5 

-125.7 

-114.4 

-87.3 

-88.0 

6 0 

-250.8 

-228.2 

-203.9 

-193.4 

-192.4 

-207.8 

-188.3 

6 2 

63.0 

77.9 

73.8 

71.1 

75.2 

93.4 

80.4 

6 4 

142.4 

114.9 

102.4 

94.6 

88.5 

68.3 

67.1 

6 6 

147.5 

108.8 

95.2 

86.2 

75.8 

54.1 

53.9 

Ce:LuAG Nd:LuAG Sm:LuAG Gd:LuAG Dy:LuAG Er:LuAG Yb:LuAG 

2 0 

-36.3 

22.7 

33.8 

43.7 

52.2 

80.9 

77.1 

2 2 

-33.0 

-20.4 

-14.2 

-9.1 

-3.5 

18.8 

17.8 

4 0 

4.6 

1.9 

-2.1 

-2.0 

-1.8 

-28.9 

-27.1 

4 2 

320.4 

255.8 

209.1 

184.5 

165.6 

188.6 

164.5 

4 4 

-207.8 

-165.1 

-132.0 

-114.5 

-105.3 

-97.1 

-83.6 

6 0 

-252.2 

-229.1 

-197.8 

-178.8 

-169.6 

-210.4 

-185.5 

6 2 

67.6 

81.3 

73.9 

68.8 

68.3 

79.1 

69.0 

6 4 

141.7 

118.7 

97.0 

84.7 

78.5 

71.9 

60.9 

6 6 

148.8 

109.0 

89.6 

75.9 

66.2 

57.8 

50.5 


V. DISCUSSION 

As seen in Fig. 1, the hybridization parameter A for R imparities in YAG and LnAG comes ont 
smaller than in manganites and LaFs. This represents a serions obstacle when trying to get the best 
GFP, as for A smaller than ~ 4 eV the calcnlation does not always yield reliable resnlts. We traced 
the problem to the wannierQO package. Galcnlation with wannierQO provides the maximally localized 
Wannier fnnctions, bnt it does not gnarantee that they will be centered on the crystal site of the R 
imparity. Indeed, for A ~ 4 eV the fnnctions are displaced for heavier R and for still smaller A the 
displacement appears for all R. As a conseqnence the symmetry is lost, all GFP are nonzero and 
for q ^ 0 they are complex. A possible remedy may be to nse recently proposed scheme by Saknma 


9 








FIG. 6: Nd:YAG. Splitting of the lowest ten multiplets by the crystal field. The solid line is only to guide the 
eye. The lowest energy of the j-ih. multiplet (j=l,...,10) for calculated (experimental) results was reduced by 
0(0), 255(248), 495(489), 722(714), 1452(1418), 1567(1533), 1694(1656), 1851(1813), 1999(1952), 2020(2088) 
meV, respectively. Experimental data were taken from Ref. [21. 


32| {symmetry adapted Wannier functions) or simpler, but less sophisticated method of selectively 


localized Wannier functions 

Despite the problem with A the agreement between calculated and experimental multiplet splitting 
is very good as shown in Figs. IM7II51 We found similar agreement also for other garnet systems. 

In the past the CFP of the R impurities in garnets were determined several times using either 
semiempirical methods or the least squares £t to the optical data. Prominent groups that adopted this 
approach are those of Gruber and Burdick, which obtained rich optical data and carefully analyzed 
them using the least squares £t. In Table IIVI we compare our CFP with those obtained by these 
groups for Nd:YAG [aij, Sm:YAG pj and Er:YAG [23 1. 

When considering one should be aware of the fact that the spectroscopic methods do not 
differentiate between the six crystallographically equivalent R sites which, however, have different 
local coordinate system. The local coordinate systems are connected by symmetry operations of the 
laid space group and to c omp are different sets of CFP corresponding symmetry operation has to be 


applied (see e.g. Refs. 


23 


Ml- 


In Table HVj are also compared the quantities Sk, which were introduced by Leavitt 3^ and which 
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FIG. 7: Sm:YAG. Splitting of the lowest ten multiplets by the crystal field. The solid line is only to guide the 
eye. The lowest energy of the j-th multiplet (j=l,...,10) was reduced for calculated (experimental) results 
by 0(0), 129(126), 284(279), 446(440), 609(603), 768(7^, 909(813), 1017(831), 1164(833), 1331(848) meV, 


respectively. Experimental data were taken from Ref. 
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TABLE II: g tensor components of the ground Kramers doublets along principal axes in YAGrR^'*' calculated 
for hybridization parameter A=5.4 eV, except for Er where A=4.8 eV. Experimental data refer to the work 
listed in the last column. Difference in per cent between calculated and experimental values is reported. 


R 


9a 



gb 



9c 


Ref. 


calc. 

exp. 

diff. 

calc. 

exp. 

diff. 

calc. 

exp. 

diff. 


Ge 

1.06 

0.91 

16 

2.54 

1.87 

36 

2.32 

2.74 

15 

m 

Nd 

2.16 

1.74 

24 

1.74 

1.16 

50 

3.52 

3.91 

10 

[Ml 

Sm 

0.59 

— 

— 

0.13 

— 

— 

0.02 

— 

— 

— 

Gd 

1.82 

1.99 

9 

1.82 

1.99 

9 

2.43 

1.99 

22 

[28J 

Dy 

0.07 

0.40 

83 

0.30 

0.73 

59 

18.9 

18.2 

4 

[29J 

Er 

8.41 

7.75 

9 

3.77 

3.71 

2 

7.50 

7.35 

2 

[29J 

Yb 

3.97 

3.78 

5 

3.83 

3.87 

1 

2.52 

2.47 

2 

[25] 
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FIG. 8: Er:YAG. Splitting of the lowest ten multiplets by the crystal field. The solid line is only to guide the 
eye. The lowest energy of the j-th multiplet (j=l,...,10) for calculated (experimental) results was reduced 
by 0(0), 812(812), 1271(1272), 1525(1525), 1895(1896), 2280(2281), 2367(2370), 2543(2543), 2755(2755), 
2801(2803) G 



FIG. 9: Ce:YAG, LuAG. Principal components of g as a function of the 
imental data (enlarged symbols on the left side) are taken from Refs. [26 


l ybr idization parameter A. Exper- 
27. 
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FIG. 10: Nd:YAG, LuAG. Principal components of 5 as a function of t 
Experimental data (enlarged symbols on the left side) are taken from Refs. 


le hybridization parameter A. 
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FIG. 11: Sm:YAG, LuAG. Principal components of g as a function of the hybridization parameter A. 
Experimental data are not available. 
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FIG. 12 : Gd:YAG, LuAG. Principal components of 5 as a function of the hybridization parameter A. 

TABLE III: g tensor components of the ground Kramers doublets along principal axes in LuAG:R^"*“ cal¬ 
culated for hybridization parameter A=5.4 eV, except for Er where A=4.8 eV. Experimental data refer to 
the work listed in the last column. Difference in per cent between calculated and experimental values is 
reported. 


R 


9a 



9b 



9c 


Ref. 


calc. 

exp. 

diff. 

calc. 

exp. 

diff. 

calc. 

exp. 

diff 


Ge 

1.23 

0.92 

13 

2.33 

1.87 

25 

2.25 

2.61 

14 

[27] 

Nd 

2.20 

1.79 

23 

1.79 

1.24 

44 

3.45 

3.83 

10 

m 

Sm 

0.65 

— 

— 

0.14 

— 

— 

0.002 

— 

— 

— 

Gd 

1.88 

— 

— 

1.81 

— 

— 

2.37 

— 

— 

— 

Dy 

0.55 

2.29 

76 

0.11 

0.91 

88 

18.74 

16.6 

13 

[29] 

Er 

6.87 

4.12 

67 

4.01 

8.43 

52 

8.66 

6.93 

25 

[29] 

Yb 

3.83 

3.82 

0.3 

4.02 

3.72 

8 

2.48 

2.57 

4 

m 


are invariant with respect to the rotation of the coordinate system: 


Sk = 


2k+ 1 


E I' 


q=—k 


1/2 


(4) 


We now turn to the magnetism, which is more susceptible to the values of CFP compared to the 
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FIG. 13: Dy:YAG, LuAG. Principal components of 5 as a function of the hybridization parameter 
Experimental data (enlarged symbols on the left side) are taken from Ref. 


A. 



FIG. 14: Er:YAG, LuAG. Principal components of 5 as a function of the hybridization parameter A. 


Experimental data (enlarged symbols on the left side) are taken from Ref. 


29. 
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FIG. 15: Yb:YAG, LuAG. Principal components of 5 as a function of 
Experimental data (enlarged symbols on the left side) are taken from Refs. 

energies in the zero magnetic field. Despite this sensitivity and the above mentioned problem with A 
smaller than 4 eV, the calculated g tensors qualitatively reflect the experimental data. The method 
will be thus particularly useful for magnetic, nontransparent rare-earth compounds, in which the 
experimental data do not allow determination of the crystal held. 

Our approach is relatively simple, it may be used by non specialists and corresponding programs, 
as well as the test example, are available on the WIEN2k web site (www.wien2k.at). 

VI. CONCLUSIONS 

The crystal held parameters and the g tensors of the ground state were calculated for all seven 
rare-earth Kramers ions substituted for Y (Lu) in YAG (LuAG). Very good agreement with the 
spectroscopic data and qualitative agreement with experimental g tensors was found. A comparison 
of calculated GFP with their counterparts obtained by the least squares ht to the optical data is fair 
and it shows that theory can help to avoid the ambiguity inherent to the least squares approach. 


he hybridization parameter A 
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TABLE IV: Comparison of CFP and the rotational invariants Sk obtained in this paper and CFP obtained 
by Gruber et al. and Burdick et al. (these parameters were rotated in the coordinate system used in this 
paper as described in the text). CFP and are in meV, hybridization parameter equals to 4.76 eV for Nd, 
6.8 eV for Sm and Er. 


Nd Sm Er 


k q this paper Ref. 21 this paper Ref. 22 this paper Ref. 23 


2 

0 

32.8 

52.3 

37.0 

54.3 

113.1 

42.3 

2 

2 

-28.7 

-19.7 

-23.8 

-13.3 

15.4 

-27.6 

4 

0 

-5.8 

-23.5 

-1.1 

-16.1 

-22.9 

-21.4 

4 

2 

253.9 

273.4 

223.2 

238.7 

183.8 

185.5 

4 

4 

-161.9 

-119.5 

-139.5 

-60.7 

-87.3 

-52.1 

6 

0 

-228.2 

-241.1 

-203.9 

-197.0 

-207.8 

-146.0 

6 

2 

77.9 

73.5 

73.8 

85.2 

93.4 

40.0 

6 

4 

114.9 

103.5 

102.4 

117.2 

68.3 

65.6 

6 

6 

108.8 

72.2 

95.2 

85.5 

54.1 

64.7 

2 

23.3 

26.5 

22.4 

25.7 

51.5 

25.7 


Sk 4 142.0 

140.9 

124.1 

116.2 

96.2 

91.1 


6 

93.8 

88.0 

83.9 

85.7 

76.4 

56.5 
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